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We present a class of exact solutions of Einstein's gravitational field equations describing spher- 
ically symmetric and static anisotropic stellar type configurations. The solutions are obtained by 
assuming a particular form of the anisotropy factor. The energy density and both radial and tan- 
gential pressures are finite and positive inside the anisotropic star. Numerical results show that the 
" " J" . basic physical parameters (mass and radius) of the model can describe realistic astrophysical objects 

' like neutron stars. 
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00 ' Relativistic stellar models have been studied ever since the first solution of Einstein's field equation for the interior 
of a compact object in hydrostatic equilibrium was obtained by Schwarzschild in 1916. The search for the exact 
solutions describing static isotropic and anisotropic stellar type configurations has continuously attracted the interest 
of physicists. The study of general relativistic compact objects is of fundamental importance for astrophysics. After 
the discovery of pulsars and explanation of their properties by assuming them to be rotating neutron stars, the 
theoretical investigation of superdense stars has been done using both numerical and analytical methods and the 
parameters of neutron stars have been worked out by general relativistic treatment. 

There are very few exact interior solutions (both isotropic and anisotropic) of the gravitational field equations 
~^ . satisfying the required general physical conditions inside the star. From 127 published solutions analyzed in Delgaty 
^ ' and Lake (1998) only 16 satisfy all the conditions. But the study of the interior of general relativistic stars via finding 
, exact solutions of the field equations is still an active field of research. Hence recently an algorithm was proposed 
O^' by Fodor (2000) to generate any number of physically realistic pressure and density profiles for spherical perfect 
^ [ fluid isotropic distributions without evaluating integrals. The gravitational field equations for static stellar models 
with a linear barotropic equation of state Nilsson and Uggla (2001b) and with a polytropic equation of state Nilsson 
^ ' and Uggla (2001a) p — kp^^^^" were recast respectively into two complementary 3-dimensional regular system of 
ordinary differential equation on compact state space, these systems were analyzed numerically and qualitatively, 
using the theory of dynamical system. Schmidt and Homann (2000) discussed numerical solutions of Einstein's field 
equation describing static spherically symmetric conglomerations of a photon star with an equation of state p = 3p. 
Recently upper limits for the mass-radius ratio are derived for compact general relativistic objects in the presence of 
a cosmological constant (Mak, Dobson and Harko 2000) and in the presence of a charge distribution (Mak, Dobson 
and Harko 2001). 

Since the pioneering work of Bowers and Liang (1974) there is an extensive literature devoted to the study of 
anisotropic spherically symmetric static general relativistic configurations. The theoretical investigations of Ruderman 
(1972) about more realistic stellar models show that the nuclear matter may be anisotropic at least in certain very 
high density ranges {p > lO^^g/cm'^), where the nuclear interactions must be treated relativistically. According to 
these views in such massive stellar objects the radial pressure may not be equal to the tangential one. No celestial 
body is composed of purely perfect fluid. Anisotropy in fluid pressure could be introduced by the existence of a solid 
core or by the presence of type 3 A superfluid (Kippenhahm and Weigert 1990), different kinds of phase transitions 
(Sokolov 1980), pion condensation (Sawyer 1972) or by other physical phenomena. On the scale of galaxies, (Binney 
and Trcmainc 1987) have considered anisotropics in spherical galaxies, from a purely Newtonian point of view. 
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The mixture of two gases (e.g., monatomic hydrogen, or ionized hydrogen and electrons) can formally be also de- 
scribed as an anisotropic fluid (Lctelicr 1980, Bayin 1982). Bowers and Liang (1974) have investigated the possible 
importance of locally anisotropic equations of state for relativistic fluid spheres by generalizing the equations of hydro- 
static equilibrium to include the effects of local anisotropy. Their study shows that anisotropy may have non- negligible 
effects on such parameters as maximum equilibrium mass and surface redshift. Heintzmann and Hillebrandt (1975) 
studied fully relativistic, anisotropic neutron star models at high densities by means of several simple assumptions 
and have shown that for arbitrary large anisotropy there is no limiting mass for neutron stars, but the maximum mass 
of a neutron star still lies beyond 3 — AMq. Hillebrandt and Steinmetz (1976) considered the problem of stability 
of fully relativistic anisotropic neutron star models. They derived the differential equation for radial pulsations and 
showed that there exists a static stability criterion similar to the one obtained for isotropic models. Anisotropic 
fluid sphere configurations have been analyzed, using various Ansatze, in Bayin (1982), Cosenza, Herrera, Esculpi 
and Wittcn (1981) and Harko and Mak (2000). For static spheres in which the tangential pressure differs from the 
radial one, Bondi (1992) has studied the link between the surface value of the potential and the highest occurring 
ratio of the pressure tensor to the local density. Chan, Herrera and Santos (1993) studied in detail the role played by 
the local pressure anisotropy in the onset of instabilities and they showed that small anisotropies might in principle 
drastically change the stability of the system. Herrera and Santos (1995) have extended the Jeans instability criterion 
in Newtonian gravity to systems with anisotropic pressures. Recent reviews on isotropic and anisotropic fluid spheres 
can be found in Delgaty and Lake (1998) and Herrera and Santos (1997). 

Very recently an analysis based on the Weyl tensor of the Lemaitre-Schwarzschild problem of finding the equilibrium 
conditions of an anisotropically sustained spherical body under its relativistic gravitational field has been revisited 
by Fuzfa, Gerard and Lambert (2001). Hernandez and Nunez (2001) presented a general method for obtaining static 
anisotropic spherically symmetric solutions satisfying a nonlocal equation of state from known density profiles, by 
assuming the condition of a vanishing Weyl tensor. This condition can be integrated in the spherically symmetric case. 
Then, the resulting expression is used by Herrera, Frisco, Ospino and Fuenmayor (2001) to find new, conformally 
flat, interior solutions to Einstein equations for locally anisotropic fluid. Dev and Gleiser (2000) presented several 
exact solutions for anisotropic stars of constant density. Some classes of exact solutions of Einstein field equations 
describing spherically and static anisotropic stellar type configurations were presented by Harko and Mak (2002) and 
Mak, Dobson and Harko (2001). The surface redshift Zg of anisotropic realistic stars has been investigated by Ivanov 
(2002). Zg cannot exceed the values 3.842 or 5.211, when the tangential pressure satisfies the strong or dominant 
energy conditions, respectively. 

In the present paper we consider a class of exact solutions of the gravitational field equations for an anisotropic 
fiuid sphere, corresponding to a specific choice of the anisotropy parameter. The metric functions can be represented 
in a closed form in terms of elementary functions. All the physical parameters like the energy density, pressure and 
metric tensor components arc regular inside the anisotropic star, with the speed of sound less than the speed of 
light. Therefore this solution can give a satisfactory description of realistic astrophysical compact objects like neutron 
stars. Some explicit numerical models of relativistic anisotropic stars, with a possible astrophysical relevance, are also 
presented. 

This paper is organized as follows. In Section II we present an exact class of solutions for an anisotropic fiuid sphere. 
In Section III we present neutron star models with possible astrophysical applications. The results are summarized 
and discussed in Section IV. 



II. NON-SINGULAR MODELS FOR ANISOTROPIC STARS 



In standard coordinates a;' = {t, r, x, (f>), the general line element for a static spherically symmetric space-time takes 
the form 

ds"" = A^{r)dt^ - V-\r)dr^ - (^x' + sin^ x#') . (2.1) 
Einstein's gravitational field equations are (where natural units 8nG = c = 1 have been used throughout): 

Ri - IrS^ = Tt- (2.2) 

For an anisotropic spherically symmetric matter distribution the components of the energy-momentum tensor are 
of the form 

Tt = {p + P±) Uiu'^ - p^S'y + {pr - p±) Xix\ (2.3) 
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where is the four- velocity Au^ = (5q, is the unit spacchkc vector in the radial direction = ^/V6l, p is the 
energy density, Pr is the pressure in the direction of x' (normal pressure) and p±^ is the pressure orthogonal to Xi 
(transversal pressure). We assume Pr 7^ P±- The case pr — p± corresponds to the isotropic fluid sphere. A = p± — Pr 
is a measure of the anisotropy and is called the anisotropy factor Herrcra and Ponce de Leon (1985). 

A term appears in the conservation equations T^.^ = 0, (where a semicolon ; denotes the covariant derivative 

with respect to the metric), representing a force that is due to the anisotropic nature of the fluid. This force is directed 

outward when p±^ > pr and inward when p± < p^. The existence of a repulsive force (in the case p±^ > pr) allows 
the construction of more compact objects when using anisotropic fluid than when using isotropic fluid (Gokhroo and 
Mehra 1994). 

For the metric (2.1) the gravitational field equations (2.2) become 



VM'4 + iV^v(4^-^-lU2(A-JL|, (.5) 



A r I \ A rA I \ 



where 1 = 4-- 

dr 



To simplify calculations it is convenient to introduce the following substitutions: 

m(r) 1 



V = l- 2x7?, x = r\'n(r) = m(r) = - / (0 d^- (2-6) 

1" ^ Jo 

m{r) represents the total mass content of the distribution within the fluid sphere of radius r. Hence, we can express 
(2.5) in the equivalent form 

For any physically acceptable stellar models, we require the condition that the energy density is positive and finite 
at all points inside the fluid spheres. In order to have a monotonic decreasing energy density p — (r^r"*) inside 
the star we chose the function r] in the form 

where ip = cqx, and oq, cq are non-negative constant. (2.8) represents an ansatz for the mass function, which has 
been previously used in the study of the isotropic fluid spheres by Matese and Whitman (1980) and Finch and Skea 
(1989). 

We also introduce a new variable A by means of the transformation 

^^(ao-co)(l+V>) ^2.9) 



We also chose the anisotropy parameter as 



[ip+l) 



where ao/co = or cquivalcntly Aq = (1 — K) /cq + K and K a constant. In the following we assume Aq > 0, 

with Aq — corresponding to the isotropic limit. Hence Aq can be considered, in the present model, as a measure 
of the anisotropy of the pressure distribution inside the fluid sphere. At the center of the fluid sphere the anisotropy 
vanishes, A(0) = 0. For small r the anisotropy parameter increases and, after reaching a maximum in the interior of 
the star, it becomes a decreasing function of the radial distance. 

Therefore with this choice of A, (2.7) becomes a hypergeometric equation, 

, , d^A IdA K , , , 
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by 



By introducing the substitution A = 1 — X"^, (2.11) reduces to a standard second order differential equation given 



(FA „ dA 



By means of the pair of transformations (Gupta and Jasim 2000) 



(2.12) reads 



§=G,X = sin., 



(2.12) 



(2.13) 



(2.14) 



where w'^ = l-\- K. 

Wc shall not present here a simple linear solution for G corresponding to uP = 0. For nonzero and with the use 
of the general solution of (2.14), (2.13) becomes 



dA 

— — = G = pi cos (w sin~^ X) + q\ sin (uj sin~^ X) , 
dX 



(2.15) 



where pi and qi are arbitrary constants of integration. 

In view of (2.12) and (2.15) and after some re-adjustments of the arbitrary constants, we obtain the general solution 
of (2.15) in the form 



A = ai 



1 — X'^ cos (u) sin ^ X) + X sin (a; sin ^ X) 
u)\/l — X^sin {uj sin~^ X^ — X cos (w sin~^ X) 



+ 



(2.16) 



where a\ and /3i are arbitrary constants of integration. 

For simplicity, the metric functions and the physical parameters p,p±,Pr and A can be expressed in terms of 6 by 

setting sin 6* = X,soX= (^^) (1 + V) = 1 - = cos^ 9. 

Therefore the general solution of the gravitational field equations for a static anisotropic fluid sphere can be presented 
in the following form expressed in elementary functions: 



A=VC{{u/-l) sin [{lj + 1)0 + B] + {lu + 1) sin [(w - 1) 6* + B]} , 



(2.17) 



V = l 



An -if 



1-K y vi + V' 



Ao-l W.„2 



1-K 



tan^ e, 



(2.18) 



CqAo^ ^0^0 



( A^-A: ^ 



COS^ 6 - 1 



V Ao-1 ) 



cos^ 9 



(2.19) 



P = Co 



Aq- if 

1-K 



3 + V> 

(1 + V')' 



CO (Aq - 1)" 

(1 - K) (Ac - K) 



Aq-K 
Ao-1 



sec2 6' + 2sec'^6' 



(2.20) 



Pr = Co 



1- Aq 

1-K 



sec 



2K 



1- Aq 
if- An 



tan 9 



u)tan{u)9 + B) — tan^ 



(2.21) 



cqAoV' , 
(1 + ^) 



(2.22) 
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where ol\ = 2vCsmB and /3i = 2^^008 5 and B and C arc arbitrary constants. 

The physical quantities ^ and describing the behaviour of the speed of sound inside the static anisotropic 
fluid sphere can be computed from the resulting line element. 

In order to be physically meaningful, the interior solution for static fluid spheres of Einstein's gravitational field 
equations must satisfy some general physical requirements. The following conditions have been generally recognized 
to be crucial for anisotropic fluid spheres Herrera and Santos (1997): 

a) the density p and pressure pr should be positive inside the star; 

b) the gradients ^, ^ and ^ should be negative; 

c) inside the static configuration the speed of sound should be less than the speed of light, i.e. < ^ < 1 and 

o<^<i; 

d) a physically reasonable energy- momentum tensor has to obey the conditions p > Pr + 2p± and p + Pr + > 0; 

e) the interior metric should be joined continuously with the exterior Scliwarzscliild metric, that is A'^[a) = 1 — 2u, 
where u = M/ a, M is the mass of the sphere as measured by its external gravitational field and a is the boundary of 
the sphere; 

f) the radial pressure Pr must vanish but the tangential pressure p± may not vanish at the boundary r = a of the 
sphere. However, the radial pressure is equal to the tangential pressure at the center of the fluid sphere. 

By matching (2.18) on the boundary of the anisotropic sphere and by denoting = cqo^, we obtain 

Vanis (a) = 1 - (^^) (t^) = 1 - 2«ani«, (2.23) 



l-K J \l + ija 

where generally u denotes the mass-radius ratio of the isotropic or anisotropic star and the suffix "a" represents the 
value of the physical quantities at the vacuum boundary of the star r = a. 

For the isotropic case, that is for Aq = 0, it is easy to show that Uiso = 2(k-i) ( i+V'a ) " Therefore the mass-radius 
ratios for the anisotropic and isotropic spheres are related in the present model by 

Uanis — 'Wjso- (2.24) 

Hence, the constants C and B appearing in the solution can be evaluated from the boundary conditions. Thus, we 
obtain 

, _^ 2_ 2(l-A-)a,„„ _ 2 (1 - A-) a,,,,,,. 

C = (1 - 2uams) {{ui - 1) sin [(w + l)9a + B] + {u; + 1) sin [(w - 1) 0a + B]}-^ . (2.26) 

Using (2.21) and the boundary condition that the radial pressure Pr vanishes at the radius r = a of the fiuid sphere, 
we obtain 



tanB 



K+Aoil-2K) 



1 + 



K+Ao{l-2K) 
u>(Ao-K) 



tan^^a — tan {ujOa) 

(2.27) 



tan {(jj9a) tan^a 



By defining tan {uj6a) = [tan {u)0)]^^^, it is easy to show that 



_ / (Ao-l)(l + V>a) 



tan 9a 



(Ao - K) 



(Ao-1) (1 + ^a) 



— 1, tan {uj9a) = tan 



wtan 



(Ao - K) 



(Ao-l)(l + Va) 



(2.29) 



The mathematical consistency of (2.28) requires the condition (Aq — 1) (Aq — ii') > be satisfied. 

In order to find general constraints for the anisotropy parameter Aq, K and Uanisi we shall consider that the 
conditions po = p{0) > 0, po = p{0) > and Rq = R{0) = po — 3po > hold at the center of the fluid sphere. 
Subsequently the parameters Aq, K and Uanis should be restricted to obey the following conditions: 
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(2.30) 



(Ao - K) VI + A' + 4/v [K - 1) cot B 



>o, 



(2.31) 



6 [{K - Ao) VTTK -2K{K- 1) cot B] 
Ro = : ^ — , > 0, 



(2.32) 



Clearly, the condition of the non-negativity of the energy density is satisfied only if (Aq — K) (1 — K) > 0. 
The general behavior of the functions po [K, Aq), po {K, Aq, Uiso) and Ro{K, Aq, Mi^o) is represented in Fig.l. 
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FIG. 1. Variations of the trace Ro = po — 3po of the energy momentum tensor at the centre of the anisotropic star (solid 

curve), of the central density po (dotted curve) and of the central pressure po (dashed curve) as functions of the anisotropy 
parameter Aq, for K — 2.6 and for a given mass-radius ratio Uiso = 0.24. 



Since cos^a < 1, from (2.28) it follows that the parameters K, Aq and ipa must satisfy the condition 

(Ao-l)(l+Va) 



(Ao - K) 



< 1. 



(2.33) 



(2.33) also represents the mathematical consistency condition for equations (2.29) be well-defined. The upper limit 
of the radius a of the anisotropic star can be obtained from (2.33) in the form: 



a < 



I 1-K 
Co (Ao - 1) ■ 



(2.34) 



(2.34) requires for the parameters K and Ao to satisfy the condition (1 — K) {Aq — 1) > 0. 

Since from (2.23) ipa can be defined in terms of the total mass Manis as ipa = Ao-K-2{i-K)M^is/a ' (2-33) gives 
the minimum allowed mass of the anisotropic star in the form 



Ao-1 
2(1 -K) 



Co a 



(2.35) 



The conditions of the non-negativity of the energy density and maximum radius of the star give the allowed 
ranges of the; parameters Aq and K. There are two ranges for these parameters making the physical requirements 
mathematically consistent, and they are given by: (i) K < 1 and Aq > 1 > K and (ii) K > 1 and Ao < 1 < For 
values of Aq and K which do not belong to these sets, the field equations do not have a physical solution. 
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III. ASTROPHYSICAL APPLICATIONS 



When the thermonuclear sources of energy in its interior are exhausted, a spherical star begins to collapse under the 
influence of gravitational interaction of its matter content. The mass energy continues to increase and the star ends 
up as a compact relativistic cosmic object such as neutron star, strange star or black hole. Important observational 
quantities for such objects are the surface redshift, the central redshift and the mass and radius of the star. 

For a relativistic anisotropic star described by the solution presented in the previous Section the surface redshift Zg 
is given by 

= (1 - 2Uanis)~^'^ - 1. (3.1) 

The surface redshift is decreasing with increasing Aq. Hence, at least in principle, the study of redshift of light 
emitted at the surface of compact objects can lead to the possibility of observational detection of anisotropics in the 
internal pressure distribution of relativistic stars. 

The central redshift Zc is of the form 

z^ = C-^l^ {(w - 1) sin [(w + 1) ^0 + -B] + (w - 1) sin [(w + 1) + -B]}"'/^ - 1, (3.2) 



where = cos ^ ^ Aq-k ■ variation of the central redshift of the neutron star against the anisotropy parameter 
is represented in Fig. 2. 




0.1 0.2 0.3 0.4 0.5 0.6 

Ao 

FIG. 2. Variation of the central redshift Zc as a function of the anisotropy parameter Ao for the static anisotropic fluid 
sphere with K = 2.2 and a mass-radius ratio Uiso = 0.24. 

Clearly, in our model the anisotropy introduced in the pressure gives rise to a decrease in tlic; central redshift. 
Hence, as functions of the anisotropy, the central and surface rcdshifts have the same behavior. 

The stellar model presented here can be used to describe the interior of the realistic neutron star. Taking the 
surface density of the star as = 2 x lO^^g/cm^ and with the use of (2.25) we obtain 

2 _ {Ao-K)yJa (3 + V>a) 
[1 - K) {1 + iPa) 

or 

23.144 



ttN = V(Ao - K) Uiso [4 {K - 1) Ui,o - 3K]km, (3.4) 

where cat is the radius of the neutron star corresponding to a specific surface density. 
For the mass M;v and anisotropy parameter A;v we find 

15.67 



Mn = {K - Ao) UisoVi^o - K) Uiso [4 {K - 1) Uiso - ^K]Mq, (3.5) 
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= 3.59502 x 1035- j,^,, ^r^" (3-6) 

In the equations above we have expressed, for the sake of simphcity, all the quantities in international units, 
instead of natural units, by means of the transformations Stt n — + pc? and A)v — > ^^Aw, where 

G = 6.6732 X IQ-^dyne cm^ 5-2 and c = 2.997925 x IQi^cm s"!. 

The variation of the anisotropy parameter Aat of the neutron star as a function of the radius aN and mass Mn is 
represented in Fig. 3. 




The quantities ( ^ ) , ( ) , ( ^ ) and ( ) are represented against the anisotropy parameter 
Aq in Fig. 4. 
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FIG. 4. Variations of the radial and tangential speeds of sound dp/dp at the vacuum boundary and at the center of the 
anisotropic fluid sphere as a function of the anisotropy parameter Ao for Uiso = 0.2 and K = 2.2: {dpr/dp)^^Q (solid curve), 
{dp±/dp)^^Q (dotted curve), (dpr/dp)^^^^ (dashed curve) and {dp±/dp)^^^^ (long dashed curve). 

The plots indicate that the necessary and sufficient criterion for the adiabatic speed of sound to be less than 
the speed of light is satisfied by our solution. However, Caporaso and Brcchcr (1979) claimed that dp/ dp docs not 
represent the signal speed. If therefore this speed exceeds the speed of light, this does not necessary mean that the 
fluid is non-causal. But this argument is quite controversial and not all authors accept it Glass (1983). 

IV. DISCUSSIONS AND FINAL REMARKS 

In the present paper we have presented the general solution of the gravitational field equations for an anisotropic 
static matter distribution. To derive the solution we have used two basic assumptions: the functional form of the 
mass function ry(r) and a specific mathematical representation of the anisotropy parameter A. In our model the 
anisotropy vanishes at the center of the star, reaching its maximum near the vacuum boundary. This behavior is 
similar, from a phenomenological point of view, with the behavior of the anisotropy parameter of boson stars (Dev 
& Gleiser 2001). Boson stars are gravitationally bound macroscopic quantum states made of scalar bosons (Jetzer 
1992). They differ from the fermionic stars in that they are only prevented from collapsing gravitationally by the 
Heisenberg uncertainty principle. Boson stars, described by non-interacting, massive scalar fields (j) with non-zero 

spatial gradients arc anisotropic objects, with = — pj. = — . Hence our A has the same qualitative 

properties as the anisotropy parameter of boson stars. 

An other possible source of the anisotropy as described by (2.10) could be an anisotropic velocity distribution of 
the particles inside the star. Most of the models of stellar structure assume that the equation of state is due to a 
single species of ideal noninteracting fermions, with an isotropic distribution of momenta. If, due to the presence of a 
magnetic field, convection, turbulence etc. the local velocity distribution is anisotropic, then the pressure distribution 
is also anisotropic. To evaluate this effect let us consider a spherical configuration of collisionless particles, which 
exhibit no net velocity flow. The radial, azimuthal and polar velocities are w^, and vg, respectively. Then one can 
show that for a system in Newtonian equilibrium A = p_L — = p [(wg) — (^r )] : where the brackets denote averages 
over the particle distribution (Herrera & Santos 1997). Moreover, {vg) — (u^) = — (Herrera & Santos 1997), 
where $ is a parameter measuring the anisotropy of velocity distributions, which for simplicity we assume to be a 
constant. Therefore 



A = p [{vl) - {vl)] = -co 



Ao 



K 



K 



3 + ^ 

(1 + V')' 



(4.1) 



But with the use of the virial theorem (the kinetic energy of particles is equal to minus half of their potential energy) 
we obtain (yl) ~ —m/r, where m is the mass of the star. Approximating the density by a constant (at least locally 
this is generally a good approximation for general relativistic compact objects) it follows that (u^) ~ — ^ —ij). 
Substituting this expression of the average of the square of the radial velocity in (3.1), by chosing an appropriate value 
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for the constant $, and by neglecting in the numerator the second order term proportional to V'^ {ip ^ r"^ /a^ << 1 
for r < a), we obtain the exact form of the anisotropy ansatz (2.10). 

The scalar invariants of the Riemann tensor are important since they allow a manifestly coordinate invariant 
characterization of certain geometrical properties of geometries, like for example curvature singularities, Petrov type 
of the Weyl tensor etc. Two scalars which have been extensively used in the physical literature are the Kretschmann 
scalars RiemSq = RijuR'"''^ = (2A + p + - + 2 (p, + + 2 {p - +A {^f and RicciSq = RijW^ = 

+ + 2A (A + 2pr), where Rij^i is the Riemann curvature tensor. 

For space-times which are the product of two 2-dimensional spaces, one Lorentzian and one Riemannian, subject 
to a separability condition on the function which couples the 2-spaces, Santosuosso et al. (1998) have suggested that 
the set 



C = {R,ri,r2,W2} , 



(4.2) 



form an independent set of scalar polynomial invariants, satisfying the number of dc;grees of freedom in the curvature. 
In (4.2) R = g^^g^'^Rijki is the Ricci scalar and the quantities ri, r2 and W2 are defined according to (Zakhary & 
Carminati 2001; Carminati, Zakhary & McLenaghan 2002) 



''I = 'PabAb' 



(4.3) 



r2=<t> 



ABAB^CC 



BBiCACA 



nb Qc Qa 



(4.4) 



W2 = *ABCD*BF* 



CD.j.EFAB 



(4.5) 



where 5^ = — ji?(5a is the trace-free Ricci tensor, ^i^^^^ denotes the spinor equivalent of Sab, '^abcd denotes 
the spinor equivalent of the Weyl tensor Cabcd and Cabcd denotes the complex conjugate of the self-dual Weyl tensor, 

^abcd = 2 ('^abcd - i * Cabcd)- 

In terms of the "electric" Eac = CabcdU^'u'^ and "magnetic" Hac = C*f^^^u''u'^ parts of the Weyl tensor, where is 
a timelike unit vector and C*^^^ = ^VabefC^^ is the dual tensor, the invariant W2 is given by (Santosuosso et al. 1998) 



W2 = ^ {^EtHlHl - E^ElEl) + ^ {H^HlHl - SE^E^^H^) . 



(4.6) 



For the case of a perfect fluid with anisotropic pressure distribution the set of invariants (4.2) is given by: 

R = p-3p^- 2A, (4.7) 



Q{p + Pr + -^] +{p + Pr-'2Af + 2{p + Pr + 2Af 



(4.8) 



r-2 



1 
512 



2A 

27[p + Pr + — 



ip- 



2Af + 2{p + pr + 2Af 



(4.9) 



2/3 



4to 

p + Pr + 2A- — 



9 / 2A 
P+Pr + 



^ „/ 2mV 2mV , ( 2mV 1, „ 

+ 2U, + — +2U-— +4 — --(p-3p^-2A) 



-\{P + Pr- 2Af - 1 (p-Hp^ + 2A)' 



(4.10) 



When A = and Pr = P± = P, we recover the form of the Ricci invariants r„ = c„(p -l-p)""'"^, c„ =constant, for 

class B warped product space-times with isotropic perfect fluid matter sources (Santosuosso et al. 1998). 

The variations of R, ri , r2 and of W2 at the center of the fluid sphere are represented, as a function of the anisotropy 
parameter Aq, in Fig. 5. 
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FIG. 5. Variations, at the center r = of the static anisotropic fluid sphere, of the curvature scalar R (solid curve), of the 
Ricci invariants ri (dotted curve), r2 (dashed curve) and of W2 (long dashed curve) against the anisotropy parameter Ao, for 
Uiso = 0.29 and K = 2.2. 

The scalars ri, r2 and 1^2 are finite at the center of the fluid sphere, ri and W2 monotonically decrease, while r2 
monotonically increases as the anisotropy parameter Aq increases. 

For a particular choice of the equation of state at the center of the star, 3prO = 3p_Lo = po and with a vanishing 
anisotropy parameter, Aq — 0, we obtain the result Uiso = 0.29 for a static isotropic fluid sphere. 

It is generally held that the trace T of the energy-momentum tensor must be non-negative. It is also the case that 
this trace condition is everywhere fulflUed if it is fulflUed at the center of the star Knutsen (1988). 

The purpose of the present paper is to present some exact models of static anisotropic fluid stars and to investigate 
their possible astrophysical relevance. All the solutions we have obtained are non-singular inside the anisotropic 
sphere, with finite values of the density and pressure at the center of the star. Variations of the physical parameters 
mass, radius, redshift and adiabatic speed of sound against the anisotropy parameter have been presented graphically. 
Because it satisfies all the physical conditions and requirements (a)-(f), the present model can be used to study the 
interior of the anisotropic relativistic objects. 
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